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1. Introduction 

Random walks are widely used throughout physics as a model for systems in which the state 
of the system can be viewed as evolving in a stochastic way from one time step to the next. 
Their properties have been extensively explored and the techniques to study them are well 
developed [1 1. In particular, at a large scale, random walks behave diffusively. 

The most-studied case is that of random walkers which have no memory of their past 
history. Many physical applications, however, call for a model in which the choice of possible 
directions for the walker's next step are given by probabilities which are influenced by the 
path it took prior to making that choice, so that its jumps are correlated; this is often called a 
persistent or correlated random walk f2l|. 

A walk with zero memory corresponds to the Bernoulli process of the usual uncorrected 
random walk. Walks with single-step memories are the most commonly studied cases of 
persistent random walks, where the random walker determines the direction it takes at a given 
step in terms of the direction taken on the immediately preceding step [1 J. Such walks - not 
necessarily restricted to lattices, as will be the case here - were first discussed in the context 
of Brownian motion |3| and fluid dynamics f4l|, and have since found many applications 
in the physics literature, most prominently in polymer conformation theory |5| and tracer 
diffusion in metals ||6l, but also in relation to the telegrapher's equation in the context of 
thermodynamics fT^. Previous works dealing with random walks on lattices with higher-order 
memory effects include that of Montroll |8|, with applications to models of polymers, and 
Bender and Richmond 0|. 

The state of the walker is thus specified by two variables, its location and the direction it 
took at the preceding step. The statistical properties of such persistent walks can be described 
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by simple Markov chains and have already been thoroughly investigated in the literature; see 
in particular ref. ifTOl . We will only provide a short review of results relevant to our purposes, 
with specific emphasis on the diffusive properties. 

The statistics of random walks with multi-step memory can in principle be analysed in 
terms of Markov chains, in a similar way to their single-step memory counterpart. However, 
the number of states of these chains grows exponentially with the number of steps accounted 
for This is the source of great technical difficulties, which are present already at the level of 
two-step processes. 

Of specific interest to us are random walks with two-step memory. Among the class of 
persistent walks under consideration, these are the simplest case beyond those with single- 
step memory, and are therefore relevant to problems dealing with the persistence of motion 
of tracer particles where the single-step-memory approximation breaks down. An example 
where this occurs is given in recent work by the present authors, on diffusion in a class of 
periodic billiard tables ifTTl . 

The paper is organised as follows. The general framework of walks on lattices is briefly 
reviewed in section [2] where we provide the expression of the diffusion coefficient of such 
walks in terms of the velocity auto-correlations. Successive approximation schemes for 
the computation of these auto-correlation functions are presented in sections [3] [4] and |5] 
pertaining to the number of steps of memory of the walkers, respectively 0, 1, and 2. Specific 
examples are discussed, namely the one-dimensional lattice and the two-dimensional square, 
honeycomb and triangular lattices, and their diffusion coefficients are computed. Some of 
the details of the computations presented in section|5]are deferred to appendices [Appendix A| 
and Appendix B Section [6]provides an alternative derivation of the diffusion coefficients of 
two-step memory persistent walks with special left-right symmetries. Conclusions are drawn 
in section |7] 



2. Diffusion on a lattice 



We consider the motion of independent tracer particles undergoing random walks on a regular 
lattice Their trajectories are specified by their initial position ro at time f = 0, and the 
sequence {vq, . . of the successive values vt e 1^. of their direction vectors at positions 
r, , where 1^. denotes the space of direction vectors allowed at site r, , which point to the lattice 
sites adjacent to r,. Here we consider dynamics in discrete time, so that the time sequences 
are simply assumed to be incremented by identical time steps T as the tracers move from site 
to site. In the sequel we will loosely refer to the direction vectors as velocity vectors; they are 
in fact dimensionless unit vectors. 

Examples of such motions are random walks on one- and two-dimensional lattices such 
as honeycomb, square and triangular lattices, but also include persistent random walks where 
memory effects must be accounted for, i.e. when the probability of occurrence of depends 
on the past history t>„_i , v„^2, ■ ■ ■■ 

The quantity we will be concerned with is the diffusion coefficient D of such persistent 
processes, which measures the linear growth in time of the mean-squared displacement of 
walkers. This can be written in terms of velocity auto-correlations using the Taylor-Green- 
Kubo expression: 




(2.1) 



where d denotes the dimensionality of the lattice and ^ is the lattice spacing. The 
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(dimensionless) velocity auto-correlations are computed as averages (•) over the equilibrium 
distribution /i, so that the problem reduces to computing 

{vQ-v„)= ^ VQ-v„^{{v(,,...,v„}). (2.2) 

As reviewed below, this can be easily carried out in the simple examples of random walks 
with zero- and single-step memories. The main achievement of this paper is to describe the 
computation of the velocity auto-correlations of random walks with two-step memory. All 
these cases involve factorisations of the measure /i({t;o, . . . ,■««}) by products of probability 
measures which depend on a number of velocity vectors, equal to the number of steps of 
memory of the walkers. These measures will be denoted by p throughout the paper 

The next three sections, sections |3] |4] and |5] are devoted to the computation of 
the diffusion coefficient for random walks with zero-, one- and two-step memories, 
respectively. 



3. No-Memory Approximation (NMA) 

In the simplest case, the walkers have no memory of their history as they proceed to their 
next position. This gives a Bernoulli process for the velocity trials, for which the probability 
measure factorises: 

n 

pi{{vQ,...,v„})=Y\p{vi). (3.1) 

1=0 



Given that the lattice is isotropic and that p is uniform, the velocity auto-correlation \2.2\ 
vanishes: 

{vo-v„) = 8„fl. (3.2) 
The diffusion coefficient of the random walk without memory is then given by 

£>NMA = ■ (3.3) 

4. One-Step Memory Approximation (1-SMA) 

We now assume that the velocity vectors obey a Markov process for which takes on 
different values according to the velocity at the previous step v„^i. We may then write 

n 

11{{VQ, v„}) = Y{P{v,\vi^x)p{vo) . (4.1) 

1=1 

Here, P(b|a) denotes the 1-step conditional probability that the walker moves in a direction 
b,given that it had direction a at the previous step. 

We denote by z the coordination number of the lattice, i.e. the number of neighbouring 
sites accessible from each site, and we denote by R the rotation operation which takes a vector 
V through all the lattice directions v,Rv,..., R^^^v. In general, the set of allowed orientations 
of V depends on the lattice site, such as in the two-dimensional honeycomb lattice. We denote 
by T the symmetry operator that maps a cell to its neighbors, which corresponds simply to 
the identity for square lattices and to a reflection for the honeycomb lattice. 
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The idea of our calculation is to express each velocity vector Vj^ in terms of the first one, 
vq, as Dj. = R'*^T^t;o> where 4 hes between and z — 1. Substituting this into the expression 



for the velocity auto-correlation {vq ■ v„), equation (2.2 1, we obtain 



VQ-v„Y\p{vi\vi_i)p{vQ)= Y, 

1=1 in i« = 



•m,- 



)Pio 



(4.2) 



vo,...,v„ 

where 

m,„,„_, =P(R'"T"-uo|R'"-'T"-i-i;o) (4.3) 

are the elements of the stochastic matrix M of the Markov chain associated to the persistent 
random walk, and p, = p{i) are the elements of its invariant (equilibrium) distribution, denoted 
P, evaluated with a velocity in the /th lattice direction. The invariance of P is expressed as 
Y,j rni.yP; — P/- These notations will be used throughout this article. 

The terms involving M in (4.2 1 constitute the matrix product of n copies of M. 
Furthermore, since the invariant distribution is uniform over the z possible lattice directions, 
we can choose an arbitrary direction for t;o, and hence write 



(t>o • Vn) = vq ■ T"t;om|" + vq ■ RT"r>om 



(«) 

'2,1 ^ 



where m|"' denote the elements of M". 

Under special symmetry assumptions to be discussed in the examples below, one has 



{VQ-V„) = (cos0)", 



(4.5) 



where (cos 0) denotes the average angle between two successive velocity vectors. It is then 
a general, well-known, result for such symmetric persistent random walks with single-step 
memory ^ that their diffusion coefficients have the form 

1 + (cos d) , , 

£>1SMA = £)nma:; 7^ (symmetric walks). (4.6) 

1 — (cos Q) 

The actual value of the diffusion coefficient depends on the probabilities P(W~^v\v), 
which are parameters of the model. Specific applications of equation ( |4.6| l are given in the 
examples below, such as shown in figure[T| To simplify the notation, we denote the conditional 
probabilities of these walks by Pj, where j = 0,...,z— 1 corresponds to the relative angle 
2n i/z of the direction that the walker takes with respect to its previous step (up to a reflection 
in the case of the honeycomb lattice). These conventions are shown in figure[2] 



(a) 



(b) 



Figure 1. Examples of walks on (a) square and (b) honeycomb lattices. Note the inversion of 
the allowed directions at neighbouring sites on the honeycomb lattice. 



Diffusion coefficients for persistent random walks 



5 












(a) 



(b) 



2 



2 



2 







(d) 








(c) 



4 



5 



Figure 2. The possible directions of motion at a given step for different lattices, relative to the 
incoming direction which is shown by the arrow, are labelled from to z — 1 , corresponding to 
the angle 2%j/z that the lattice direction j makes with respect to the reference direction, up to 
a reflection inversion in the case of the honeycomb lattice, (a) One-dimensional lattice (z = 2); 
(b) honeycomb lattice (z = 3); (c) square lattice (z = 4); (d) triangular lattice (z = 6) . 



4.1. One- dimensional lattice 

The simplest case is that of a regular one-dimensional lattice. In this case, each site is 
equivalent, and so T is the identity. Each velocity vector Vj^ has only two possible values, 
V and —V, so that R is a reflection. We denote by Pq = ^(I'k) the probability that the random 
walker continues in the same direction at the next step, and by Pi = P(— v|v) the probability 
that it reverses direction. 

It is easy to check that the velocity auto-correlation \A.A\ yields \A.5\ : 



4.2. Two-dimensional square lattice 

On a two-dimensional square lattice, Vk can take four possible values, and each lattice site is 
again equivalent. Thus T is the identity and R can be taken as an anticlockwise rotation by 
angle n/2. We denote by Pq = P{v\v) the probability that the particle proceeds in the same 
direction as on its previous step, by Pi =P{Rv\v) the probability that the particle turns to the 
left relative to its previous direction, by P2 = P{Vi?v\v) the probability that it turns around, 
and by P3 = P{?l?v\v) the probability that it turns right. 





(4.7) 



(4.8) 
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From (4.4 1, the velocity auto-correlation (t>o • is given by 



1,1 



m 



(") 



1,3 ■ 



(4.9) 



(4.10) 



The transition matrix M given by (4.3 1 is thus the following cyclic matrix: 
/ Pq Pi P2 Pi \ 

Pi Pi Po Pi ■ 

V Pi P2 Pi Po ) 

in) 

To calculate the elements m^- y of the powers M", it is possible to compute the eigenvalues 
and eigenvectors of M and then decompose it as M = Q • L • ' , where L is the diagonal matrix 
of the eigenvalues of M and Q the matrix of its eigenvectors. This procedure is, however, not 

(n) J I 

necessary here, since we only require the combination of the m; which appears in (4.9 1. To 



proceeed, we label the distinct entries of M" as a-" , using the fact that M" is also cyclic if M 
is: 

/ An) (n) 



Writing M" = 



M" 



/in-l 



An) in) (n) (n) 
d3 d4 d[ 

(n) Jn) in) (n) 



\ a 



(4.11) 



'1 



we can exploit the particular structure of the matrix to reduce it from 



a (4 X 4)-matrix to a (2 x 2)-matrix, by considering the following differences: 



'1 



Po-Pi Pi -Pi 

Pi -Pi P0-P2 

P0-P2 Pi -Pi 

Pi -Pi P0-P2 



n-l 



Jn-l) 
^1 


Jn-l 
83 

An-l 


Jn-i) 








-Pi J 



(4.12) 



The velocity auto-correlation (4.9 1 is thus given by 

(vQ-v„) = a|"'-a^"' 



= (1 0) 



PO-P2 
Pi -Pi 



(4.13) 

1 



P0-P2 Pi -Pi 
Pi- Pi P0-P2 

Summing the previous expression over all n, and using the fact that L"=o A" (I — A) 
for a matrix A, where I is the identity matrix, we obtain the diffusion coefficient (2.1 1 as 

\ / I — ^-Vi -I- ^-S h'^ — h'-t \ ^ 

1+2 1 



D 



ISMA 



= D 



NMA 



I-P0+P2 Pi -Pi 
Pi -Pi I-P0+P2 



P0-P2 
Pi -Pi 



(4.14) 



For a symmetric process in which Pi = P3, which is often imposed by a symmetry of the 

, viz. 

(4.15) 



physical system, equation (4.5 1 holds, and the diffusion coefficient takes the form (4.6 1, viz 

I+P0-P2 



D 



ISM A - 



D 



NMA 



I-P0 + P2 



If, however. Pi 7^ P3, then equation (4.6 1 is no longer valid. Instead, we have the more 
complicated expression 

l-(P0-P2)2-(Pi-P3)2 



DlSMA — Dj^MA- 



(4.16) 



-il-P0 + P2r + {Pl-PiV 

Such asymmetric walks are the lattice equivalent of the continuous-space models of persistent 
random walks with chirality considered in L12J . 
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4.3. Two-dimensional honeycomb lattice 

On the two-dimensional honeycomb lattice, shown in figure [ijb), each site has z = 3 
neighbours. Here, R is taken to be a clockwise rotatiorjjjby angle 2ti/2> and the arrangement 
of neighbours differs by a reflection T. We denote by Pq = P{—v\v) the probability that the 
particle turns around. Pi = P{—R'^v\v) that it turns left relative to its previous direction, and 
P2 = P{—Rv\v) that it turns right. The transition matrix M is thus 

/ Pa Pi Pi \ 

M= Pi Po Pi . (4.17) 
V Pi Pi Po J 
Proceeding as with the square lattice, we let 



M" 



(4.18) 



and obtain the matrix equation 

^0 Pi Pi 

Pi Pq Pi 

Pi Pi Pq 

Pq Pi Pi ' 
Pi Pq Pi 
Pi Pi Pq 

The velocity auto-correlation (|4.4|i is thus 







-aW- 


fi[a("'+ar' 




fi[aW+aW 












2 



[VQ ■ Vn) = - a 

1 

^ 2 



4")] 



1 



Pq Pi Pi 
Pi Pq Pi 
Pi Pi Pq 



n-l 



Hence the diffusion coefficient p.l| i is 

i'lSMA = D^MA 



i[ar-) 



2 

(n-l) 



[a 



(n-l) 



-a'"-')] 
-a'"-')] 

-,(n-l)l 



1+1 



1-Pq -Pi -Pi 
-P2 1-Pq -Pi 
-Pi -P2 1-Pq 




(4.19) 



(4.20) 



(4.21) 



Again, in the case of an isotropic process for which P2=Pi= Ps ('symmetric'), equation 
(4.5 I holds, and substituting P^ = (1 —Pq)/2 gives the following expression for the diffusion 



coefficient: 



D 



ISM A - 



1-Pq + Ps 



- D^MA 



3(1 -Po) 



(4.22) 



1+Pq-Ps ^ l+3Po 

For an asymmetric process for which Pi ^ P2, defining the symmetric and antisymmetric parts 
Ps = (Pi + P2)/2 and Pa = (Pi - P2) /2, we instead obtain 

l-3p2-(Po-Ps)2 



= D^MA 



3P2 + (l+Po-Ps)2' 

3(1-Po)(1 + 3Pq)-12P2 
(l+3Po)2 + 12p2 



(4.23) 



i We take a clockwise rotation as opposed to an anticlockwise one in the other examples so that the lattice directions 
are still labelled anticlockwise. 
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4.4. Two-dimensional triangular lattice 

Persistent walks on a triangular lattice were made popular by Fink and Mao ifTSl in connection 
to tie knots. Here, each site has z = 6 neighbours, so that R is an anti-clockwise rotation 
by angle k/3. Following our convention [see figure |2](d)], we denote by Pq = P{v\v) the 
probability that the particle moves forward. Pi = P{Rv\v) that it moves in the forward left 
direction relative to its previous direction, and similarly for Pj, P3, P4, and P5. 
Proceeding along the lines of the previous subsection, we let 



/ 



M" 



Jn) 

Jn) 
^1 



,{«) 



\ An) An) 
\ ^2 ^3 



(4.24) 



in) J 



in terms of which the velocity auto-correlation ( |4.4| i is given by 



(4.25) 



A computation similar to that of equation (4. 19 1 yields 



{vo-v„) = [^l ••• 0^ 

\Pi P2 ■ ■ ■ Po / 
The diffusion coefficient ( |2.1[ ) is then 



/ Po Pi ■■■ P5 \ / Po-Pi+k(Pi-Pi-P*+Pi) \ 
P5 Pq ■ ■ ■ P4 



1+1 ••• 



Pl-P4+i(P2-P3-P5+Po) 



\ P5-P2+iiP0-Pl-P3+P4) J 



( 1-Pq -Pi 
-P5 l-Po 



(4.26) 



-P4 



\ -Pi -P2 ■■■ l-Po J 



(P. 



P3 + i{Pi-P2-P4 



+P5) \ 



Pi-P^+kiPi-Pi-Ps+Po) 

V Pi-Pl+kiPo-Pl-P^+PA) ) 



(4.27) 



4.5. d-dimensional hypercubic lattice 

The case of a hypercubic lattice in arbitrary dimension d with coordination number z — 2d 
can also be treated, provided that the same probability P^ is assigned to scattering along 
all new directions which are perpendicular to the previous direction of motion. We then 
have Po+Pv/2 + 2(ii — l)Ps = 1, and the invariant distribution of velocities p{v) = 1/ {2d) is 
uniform. We then recover the expression (4. 15 1 for the diffusion coefficient in this case. 



5. Two-Step Memory Approximation (2-SMA) 

We now turn to the main contribution of this paper, namely the calculation of the diffusion 
coefficient for persistent random walks with 2-step memory, for several two-dimensional 
lattices. 
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We thus assume that the velocity vectors obey a random process for which the probabiHty 
of v„ takes on different values according to the velocities at the two previous steps, \ and 
v„_2, so that we may write 



ll{{vQ,...,v„}) =]~[P(t;,|ti,_i,'i;,_2)/5('yo,'Ui) ■ 

1=2 

The velocity auto-correlation \2.2\ function is then 



(5.1) 



(5.2) 



(t>0-f„)= ^ VQ-V„Y\P{Vi\Vi^\,Vi^2)p{vo,Vi) . 
{d„,...,do} '=2 

The calculation of these correlations proceeds in the most straightforward way by 



transposing the calculation leading to equation (2.2i to the two-step probability transitions 



as characterizing the probability transitions of a two-dimensional Markov chain. Considering 
a lattice with coordination number z, the state of the Markov chain is a normalised vector of 
dimension z^. The time evolution is specified by the (z^ x z^) stochastic matrix M with entries 

'^<=(M-i).+/2,;=(;i-i)z+i2 ^ d„,j,P{R'^JMW^Jv, W^v) , (5.3) 

where ii, ij, ii, and i2 take values between 1 and z- Denoting by P the invariant distribution of 

this Markov chain, i.e. the z -dimensional vector with components p,- such that L}=i '^ijPj — 
Pi, equation (|5.2|i becomes 



{vo-v„) 



'0.'" = 1 



(") 



m; 



where j'"' is the (z^ 



X z 



(o---i„,io P'o ! 

'0.'" = 1 

matrix with elements 



j(n) 



e •T"e- 



(5.4) 



(5.5) 



'=('1 - 1 )z+'2 J=ijl - 1 )z+;2 

and eic denotes the unit vector along the ktii lattice direction. 

Using the symmetries of the problem and writing Pj/^ = P{R'^JR-'Tv\R^Tv,v) for the 
conditional probability of turning successively by angles 2n j/z and 2nk/z with respect to the 
current direction (with reflection by T where needed), we define = exp(2/;r/z), where / 
denotes the imaginary unit, / = a/— T, and show, through the examples below, that equation 
\5A\ reduces to the general expression 



V 



Poo 
'-^Poi 



Pm 
^Poi 



Pio 



Pio 
<l>Pn 



Pz-l.O 
<t>Pz-lA 



\ 



I 



Diag(l,0,...,0''" 



'P-^ 



n-1 



Diag(l,(/)- 



Pi 



(5.6) 

where Diag(l , . . . , (^^~^) denotes the matrix with elements listed on the main diagonal and 
elsewhere. Here, cr is a sign factor which is —1 for the honeycomb lattice and +1 for the 
other lattices, and reflects the action of T. Note that the second term in the summation is the 
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complex conjugate of the first, so that the result is real. The diffusion coefficient for persistent 
random walks with 2-step memory is therefore 

■D2SMA 



^'nma 



= 1 + CJZ 1 



1 X 



(5.7) 



-<l>Poi 
i-Poo 



-Pio 



-Pz-l.O \ 
-<l>Pz-lA 



-1 



lDiag(l,0,...,0^^-i) 



1-4 



+ 



-1 



\ -0'-=/>o.c-i -4 



'Pi,z- 



Diag(l,0 



-1 



( 

V p. 



5.7. One-dimensional lattice 

We first consider the simplest case, namely the one-dimensional lattice. The stochastic matrix 
M from equation (5.3 1 is the 4x4 matrix 

/ P(v|v,v) P(v|v,-v) 



M 



(5.8) 





P(-v|v,-v) P(-v|v,v) 

P(-v|v,v) P(-v|v,-v) 

V Piv\v,-v) Piv\v,v) J 

/ Poo Pio \ 

Pii Poi 

Pqi Pn 

V Pio Poo J 

where Poo + P01 = 1 and Pio + Pi 1 = 1. 

Considering equation (5.4i, we compute the invariant distribution of M, which is the 
vector P whose components correspond to the four states p{v,v), p{—v,v), p{v,—v) and 
p{—v,—v). Given that we must have p{—v,—v) = p{v,v) and p(v, — v) = p{—v,v), the 
equilibrium distribution is obtained as the solution of the system of equations 

piv,v) =Poo/?(v,v)+Pio/7(-v,v), 



givmg 



p{v,v) + p{-v,v) = ^, 



Pi = P4 =f(v,v) 



(5.9) 



Pio 



92 = P3 =p(-v,v) = 



2[l-Poo + Pio] ' 
1-Poo 



(5.10) 



2[l-Poo+Pio] 



The matrix J*"^ equation (5.5 1, is here the same for all n, and has the expression 
/ 1 -1 1 -1 \ 



J = 



1-11-1 
-11-11 

V -1 1 -1 1 / 



(5.11) 
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The velocity auto-correlation (|5.4|l is thus 



m 



(n-l) 
1,1 



— m 



1,4 



m 



3,4 



m 



(n-l) 
3.1 



,(«-l) 
'1.3 



m 



(n-l) 
1,2 



'3,2 



Pi 

,("-l) 
3,3 



2 11 



m 



(n-l) 
1,1 

(n-l) 



'3,4 ^"^3,1 

(n) 

Since M", the nth power of M, has the symmetries of M, its entries m. y are such that 



(n-l) ^(n-1) 
1.4 '"1,3 
(n-l) _(n-l) 



P2j , 

(n-l) 
1,2 

(n-l) 



Pi 
P2 



(5.12) 



n 

^li = 


n 


^1 ' 


(n) _ (n) ^ , {« 
'^2,1 ~ '"3,4 — "A 


(n) 
^1.2 = 


(n) 


a*"' 

^2 ' 


(n) _ (n) ^ , (n 
•"2,2 ~ '"3,3 — °3 


(n) 


(n) 


An) 
^3 , 


(n) _ (n) ^ , (n 
•"2,3 — ^"3,2 — "2 


(n) 

'^lA = 


(n) 


Jn) 
^4 ' 


(n) _ (n) ^ , (n 
'^2,4 ^ '^3,1 — "l 



(5.13) 



Writing M" = MM" , we obtain two separate sets of equations for 2 x 2 matrices, one 
involving a j""* , 34" and b^"' , b^"'' , and the other involving aj"' , ag"'' and bj"' , bg"' : 



Jn) Jn) 
b(n) 

a'"' a*"' 

J") J") 

°2 '^3 



^"00 Ao 

^"01 i^ll 

Pao Pio 

Poi Pii 



(n-l) (n-l) 

a I 

.{n-l) .(n-l) 

"4 "1 

g(n-l) g(n-l) 



(5.14) 



Note that these equations do not have a simple recursive form, since the elements of the 
matrices on the two sides do not appear in the same places. However, taking the differences 
al"' — a^"-*, 33"'' — aj"', b^"' — bj"', and bj"' — b3"-', we obtain the recursive system 



Jn) J") A") 

d J d 



An) 



_Jn) 
3 ^^2 

bl"' - bl"' 



^'oo Pio 

-Poi -Pi I 

Poo Pio 

-Poi -Pi I 

Poo Pio 

-Poi -Pi I 

Poo Pio 

-Poi -Pi I 



(n-l) (n-l) (n-l) 



'1 



— a 



b(n-l)_b(«-l) b("-')-b("-') 



— a 



(n-l) 



n-l 



n-l 



al')-a 
u(i) 



(1) Ji) Jl) 

ai a 



4 °3 

b(/' b(') 



2 

b(^' 



Poo -Pio 
-Poi Pi I 

1 

-1 



(5.15) 



Plugging this equation into equation (5.12 1, we obtain 

:2fl 1 



[vo ■ v„) 



n-l 



(5.16) 



^bo Pio y " / Pi 
-^01 -Pii ) \ -P2 
This is equation ( |5.6[ ), where = exp(2/;r/2) = —1. The diffusion coefficient is therefore 
given by equation \5.1): 



D2SUh 

,Dnma 



= 1+4 1 1 



1-Poo -Pio 
Poi 1+Ai 



Pi 

-P2 
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It is a function of the two parameters Poo and Pio; when these are equal, the process reduces 
to a walk with single-step memory, and the diffusion coefficient to that of the single-step 
memory approximation ( |4.8| l, as it should. 

5.2. Two-dimensional honeycomb lattice 



For the two-dimensional honeycomb lattice, the stochastic matrix M of equation (5.3 i is a 
(9 X 9) matrix with the following non-zero entries: 

mi,i = m5,5 = mq_q = P{v\-v,v)=PQQ, 

mi,2 = m5,6 = m<)j = P{R^v\ - R^v.v) =Pio, 

mi,3 = m5,4 = mg.g = P(R?;| - Rv,v)^P2q, 

m4,i = mg.s = mi,9 = P{Rv\-v,v)=P()2, 

m4,2 ^ mg.e = m^j = P{v\ - R^v,v) = P12 , (5.18) 
m4,3 = m8,4 = m3,8 = P(R^t)| - Rv,v) =^22, 

Poi, 



m?,! = m2,5 = m(,^) = P{R^v\-v,v) ^Pq 



m7,2 m2,6 = mej = P{Rv\ - R^t;,v) =Pii , 
m7,3 = m2,4 = m6,8 = ^'(■"1 - Rv,v) = P21 ■ 
Given the three constraints 

Pio + Pu+Pi2 = l, (5.19) 

^20 + ^^21+^22 = 1, 

the actual number of independent parameters is six. 

The invariant distribution P with components p, can be written in terms of the three 
probabiUties p{v, —v), p{v, ~Rv), p{v, —R^v), 

Pl = P5 = pg = p{v,-v), 

92 = P6 = Pi ^ p{v,-Rv) , (5.20) 

P3 = P4 = P8 = P('W,-R^'U), 

which are solutions of the system of equations: 

p{v,-v)=PoQp{v,-v)+PiQp{v,~R^v)+P2Qp{v,-Rv), (5.21) 
p{v,-R^v)=Poip{v,-v)+Pnp{v,-R^v)+P2ip{v,-Rv), (5.22) 

T 1 

p{v,-v)+p{v,-Rv) + p{v,~R^v) = -. (5.23) 
The matrix (15. 51) has the block structure 



/ Bi Bi Bi \ 

j(«) = B2 B2 B2 , (5.24) 

\ B3 B3 B3 J 

where 

-1 1/2 1/2 \ / 1/2 -1 1/2 \ / 1/2 1/2 

Bi- I -1 1/2 1/2 ,B2= 1/2 -1 1/2 ,63= 1/2 1/2 -1 | . (5.25) 
-11/21/2/ \ 1/2 -1 1/2 / \ 1/2 1/2 
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'i.i 



'1.5 



1 («) 

— mi i - 

2 ^'^ 



m 



7,9 



Pi 



^4,7 

(n) 



— TO- 

2 

1 (n) 
2-^4,4- 



7,6 " 
'4,8 



(n) 


1 (n) 


"^1.9- 


2<i 


(n) 






2 


1 («)" 




2 ^'^J 


P2 + 



(") ^ (n) ^ («) 

m;,5-2"'4,9-2'^7,i 



1 



m 

1 

- -1 
2 

1 

2' 



1.7 



2 -m 



2-4.6 

1 (n^ 



'1,3" 



'1,4' 



-m 



1,8 



P3 



m 



(n) 1 ^(n) 



3(-l) 
1 



1,1 



m 



1 («) 



7,9 

(n) 1 («) 

m^ (; 2 m^ g 



■■') 


X 






1 («) 

-2"^1,2 


1 (") 
-2"^1,6 


, (n) 
+ "^1,7 


1 (n) 
-2"^1.3 


1 («) 


, (n) 


1 («) 


(n) 


(n) ■ 


1 («) 


1 («) 

5<7 


1 ' (") 



(5.26) 



(«) 

1,4 
1 («) 

4,5 ~ 2 "'4,9 "'4,2 ~ 2 "'4,6 ~ 2 "'4,7 ^ 2 " '4.3 ~ 2 "^4,4 + "^4,8 

Proceeding along the lines of the computation presented in subsection |5.1| we obtain a 





^PiX 


1 


P2 




^ P3 / 



set of recursive matrix equations (A. 9 1 involving the coefficients of M". We refer the reader 
to Appendix A for the details of this derivation. 

We note that the coefficients which appear in equation \5.26\ satisfy identities such as, 
for instance, 

1 



in) 1 (n) 1 (n) 1 



'1,1 



1,5 



2,-7r/3, 



'1,9 



(5.27) 



Thus, letting 
( |5^ to find: 



exp(2/;r/3), we can combine the results of equation (A.9l with equation 



K.^o) = (-ir^(i 1 1) 



^00 PlQ 

(pPoi (pPn m\ 
0^Po2 fPn ^^Pii 



Poo Pw PlO 
02poi (p^Pn (p^Pii 

(j)P02 (j)Pn ^P22 



D2SMh 



1-3 




(5.28) 



This is equation ( |5.6[ l. The diffusion coefficient ( |2.1[ l is therefore given by \5.1) , which is here 

-1 



(5.29) 



Given a symmetric process for which left and right probabilities are equal, but the 
probabihty of a left-left turn is different than that of a right-left turn, we let 

2 ' (5.30) 



^02 = ^-01 ^Pb. 



PlO = P20 = Psh:Pu - ^^22 = ^^88,^^12 = ^21 = 1 ~Psb-Ps. 
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Carrying out the matrix inversions in ( |5.29| l, we find the diffusion coefficient, 

n _n 3(l-Poo)(l+PQo-fib)(2-Psb-2Pss) ,,,,, 

Z)2SMA - ^NMA^^_^^^^^^^^^j^^^^^^^^^^_g^^^^^2(l+Poo)Pss-4p2] ' ^^'^'^ 

If we further assume complete left-right symmetry and identify the probabilities of left-left 
turns and left-right turns, thus letting 

1-nb 



Pn = Pn =Pn =P22 = 



3{l-Poo){l+Poo-Psb) 



equation (5.31 1 simplifies to 

''^^^ = ''^^^l-^oo + ^sb)(l+Poo + 2nb) 
5.3. Two-dimensional square lattice 



(5.32) 



(5.33) 



For the two-dimensional square lattice, the matrix M of equation (5.3 1 is the (16 x 16) matrix 
with the following non-zero entries: 



mi,i 


= m6,6 = 


= mil. 11 = 


= mi6,i6 


= P{v\v,v) =Poo, 


mi,2 


= mej = 


= mil. 12 = 


= mi6,i3 




mi,3 


= rn6,8 = 


= mii.9 = 


mi6,i4 = 


= P{R^v\R^v,v)^P2o, 


mi,4 


= m6,5 = 


= miijo = 


= mi6.i5 


^P{Rv\Rv,v)^P30, 


ms.i 


= rnio^e 


= mi5.ii 


= m4.i6 


^P{Rv\v,v)^Po3, 


m5,2 


= mioj 


= mi5,i2 


= m4.i3 


= Piv\R\,v)^Pn, 


m5,3 


= rnio.8 


= mi5.9 : 


= m4j4 = 


^P{R\\R^V,V)=^P23, 


m5,4 


= mio.5 


= mi5,io 


= m4,i5 


= PiR^v\Rv,v)=P3i, 


mg.i 


= rni4^6 


= m3ji = 


= mg.ie = 


= P{R^v\v,v)^Po2, 


m9,2 


= mi4j 


= m3.i2 = 


= m8,i3 = 


= P{Rv\R\,v)^Pn, 


m9,3 


= mi4.g 


= m3 9 


mg.i4 = 


P{V\R^V,V)^P22, 


^9,4 


= rni4 5 


= m3.io = 


= mg.is : 


^P{R\\Rv,v)^P32, 


mi3,i 


= m2,6 = 


= mv.ii = 


mi2,i6 = 


= P{R\\v,v)^Poi, 


mi3,2 


= m2,7 = 


= m7.i2 = 


mi2,i3 = 


= P{R\\R\,v)^Pu, 


mi3,3 


= ni2,8 = 


= m7,9 = 


mi2,i4 = 


P{Rv\R^v,v)^P2u 


mi3,4 


= m2,5 = 


= mv.io = 


mi2,i5 = 


= P{v\Rv,v)=P3i. 



The matrix (15. 51) is here 



j(«) 



f Bi Bi 

B2 B2 

-Bi -Bi 

y — B2 — B2 



Bi 

B2 
-Bi 
-B2 



Bi 

B2 
-Bi 
-B2 J 



where 



B, = 



/ 1 


-1 








1 





-1 \ 


1 


-1 





, B2 = 





1 





-1 


1 


-1 








1 





-1 


V 1 


-1 


0^ 




[0 


1 





-1 / 



(5.34) 



(5.35) 
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The invariant distribution P will not be written explicitly here. Due to symmetry, only 4 
of the 16 components are distinct. These four components are most simply computed as the 
invariant vector of the following (4 x 4) matrix: 



^bi Pn Pn P31 
P02 Pn P22 P32 



normalised so that pi + p2 + P3 + P4 = 1/4. 

In terms of these quantities, the velocity auto-correlations are found to be 



\ 


( P' ^ 








P2 




P2 




P3 




P3 




V P4 / 




^ P4 / 



(5.36) 





= 4 


(n) 
^1,1 - 


in) 

ml.n ^ 


(n) 


(") 


(") , 


(n) 

S,ii - 


(n) 
-"ll3,6^ 


(n) 
'"ll3,16 


Pi 










+4 


(n) 
-"^1.7 


, (n) 


, («) 
+ "^5,2 


(«) 
- "15,12 


, (n) 
+ "19.7 


(n) 


(") 
- "Il3,2 


, (n) 
+ "113,12 


P2 






+4 


(n) 


, (n) 


(«) 


(n) 
^"15,14- 


(") 

h j - 


(n) , 


(«) 
"113,8- 


(") 
"ll3,14 


P3 






+4 


(«) 


(«) 

^1,15- 


(") , 
"15,4 + 


(") 
"15,10- 


(n) , 
"19,5 + 


(«) , 
"19,15 + 


(«) 
"l l3 ,4 - 


(n) 
"ll3,10 


P4, 



/ 


(") 

"ll.l- 


(«) 

-"ll,ll 


(n) 
"I1.I3- 


(") 
- "ll,7 


(«) 
"I1.9- 


"ll,3 




(n) 
"ll3,16 


(n) 
- "Il3 6 


(n) 
"I13.I2 


(") 
- "Il3,2 


(«) 
"I13.8 - 


(«) 
"ll3,14 




(«) 

S,ii 


(n) 

- "19,1 


(n) 
"I9.7- 


(«) 
S.13 


(«) 
"19,3- 


(«) 
"19,9 


V 


(«) 

"15,6- 


(") 
-"l5,16 


(n) 
"15,2- 


(n) 
"15,12 


(«) 
"15,14- 


{«) 
-"15,8 



4(^1 1 1 1) 



Using the results of |Appendix B| this expression reduces to 



(n) 
Hlj ■ 
,(") 
'13,4 ■ 

(n) 
9,15 
(«) 



'1,15 
'13,10 



m 



5,10 



m 



5,4 



{vo ■ v„) 



2 1111 



/ ^00 
iPoi 

-Pq2 

V -'■^'03 



Pw 
iPn 

-Pn 
~iPn 



V 



Poo 
-iPoi 
-IPq2 

iPo3 



Pio 
-iPn 
'IPn 

iPi3 



P2Q 

-iP2\ 

-IP22 

iP2i 



P2Q P30 \ 

iP2\ iPsi 

-P22 -P32 

-/P23 -'^'33 / 

n-1 



l-l 



/ 1 \ 
0/00 
0-10 
V -i J 



( PI \ 

P2 
P3 
V P4 / 



(5.37) 



^^30 

-/P31 
-li'32 
/P33 



1 














/Pi\ 





— / 












P2 








-1 









P3 











-i ) 






^ P4 / 



which is equation \5.6\ . Substituting this result into equation \2.\) , we recover the diffuson 
coefficient given by equation (5^i : 

■D2SMA 



^NMA 



= 1+4 1111 



(5.38) 



/ l-Pao 
-iPoi 
P02 
iPo3 



V 



1-^00 
'■^01 
^^02 

-iPo3 



-Pio 
l-iPu 
Pn 
iPn 

-Pw 
I + /P11 

^^12 

-iPn 



-P2Q 

-iP2\ 

I+P22 
iP23 



-P3Q 

-iPi\ 

P32 
I+/P33 / 

-^20 -^30 \ 

/P2I '^31 

l+^'22 ^^32 

-iP23 -iP33 ) 



10 





^ 






i 















-1 









\ 





/ 






/I 









/p. \ 


-/ 










P2 





-1 







P3 


\ 





i J 




I P4 / 



For a symmetric walk, we substitute 

l-Poo-^'o2 



^"03 



0)1 



^23=^21 



I-P2Q-P22 



(5.39) 
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and write 

P33 = Pll=Pss, Pw = PiQ=P,u 

Pl2 = P32 =Psb, P3I = ^3 = 1 — ^sf — Pss — Psb i 

in terms of which the diffusion coefficient is foimd to be 

DisMA = Dnma 1 1 - 4^00 + 2Poo - 3Psb + ePooPsb - SPoo^sb + Psb - 2Poo/fb + P^Pi 
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(5.40) 



+3Psf - 6P02Psf - 2PooPsf + 2P02P00Psf - Poo^sf + 2PsbAf + 3Po2P,hPsf 

+3PooPsbPsi - PoiPooPsbPsf + PooPsbPsf - 2PlPsi - 2PooPlPsf + Psf 

+P02Psf + 3PooPsf ~ ^Ol^OoPsf ~ 4PsbPsf + 2Po2PsbPsf ~ 2PooPsbPsf ~ 

+2P024 - 2 [ - (- 1 + Poo)^(-l + Psb) + [1 + Po2(-3 + Poo) - 3Poo 

+2(1 +Poo)Psb]Psf-2(-l +Po2)Psf]^ss+P|o(^02-Psb)[Po2(-2 + Psb 

+Psf + 2Pss) - 2Psb (- 1 + Psb + Psf + 2Pss)] + P20 [( 1 + P00)Psb - P02Psf] 
X [3Psb + Psf + 2Pss + Po2(-2 + Psb + Psf + 2Pss) 



-2Psb (Psb + Psf + 2Pss 



^22 



2 — Psb + 3Psf ~ 2Pss 



-PoV-2 + Psb + Psf + 2Pss) - Psf [Psb + 2PsbPsf + 2Pss 
-Psf(l -2Psf-4Pss)] +Poo[Psb(2 + 3Psf) -4(1 -Pss) 

+Psf (-2 + 3Psf + 6Pss)]] - P22 {Psi [(- 1 + Poo)' - 2(1 + Poo)Psf] 

+P02Psf[-6 + Psf + 6Pss -Poo(-2 + Psf + 2Pss) +2Psf(Psf + 2Pss)] + 



Psb 



Psf (3 + 3Po2 - 2Psf + 2Po2Psf - 4Pss) + 2(- 1 + Pss) 



+Po2o(-2 + Psf + 2Pss) -Poo[-4 + Po2Psf + 2p2f + 4(1 +Psf)Pss] 



+P20 



P02 [4 - 2Psb + 2Psf - 4Pss + 2Poo(-2 + P^b + Psf + 2Pss) 



-3Psf (Psb + Psf + 2Pss)] + Psb[-4 + Psb - 3Psf + 2Pss 

+4Psf (Psb + Psf + 2Pss) - Poo(-4 + 3Psb + 3Psf + 6Pss)] 



I [1 + Psb - P00(1 + Psb) + P20(-P02 + Psb) + P22(- 1 + PoO - Psf) 

+Psf + P02Psf] [2(- 1 + Poo)Psi - 2P02Psf + Psf [1 + P22 - Poo - P22P0O 
+P02(2 + P20 - 4Pss)] + 2[1 + P22 + P20P02 - (1 + P22)P00]Pss 
+Psb [3 + P22 + P20P02 - P22P0O - 2(1 + Po2)Psf - 4Pss 

+Poo(-3 + 2Psf + 4Pss)]]|. (5.41) 

If we further assume complete left-right symmetry, so that 

P33 = Pn = P31 = Pi3 = , (5.42) 
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then the diffusion coefficient becomes 



^2SMA = ^NMA 



-1- 

•d2 



'P22 ' 



- 2P22P20P02 ^ 



PhPoO + P^b-PllPsb 



2Pqo 

PloPsb-^PllPloPsb 



2Pl2Pm-2P22Pn)PQ2P0Q 



-P^0P02Psb-2PQ0Psb 
+ P20P02P00Psb 
-3P02Psi + P20P02Psi 



2P22P00Psb-P20P0oPsb 
-P22PooPsb-3Psf 
P22P20P02Psf 



3P22Psf 



P20P02Psb 
P22P20Po0Psb 
3P02Psf 



P20PQ2Psf- 



P00Psf-P22P00Psf 



-P02P00Psf + P22P02PQ0Psf 



l+P 



22' 



22^ 
2 r}2 



2P22P20P02+Pi0P02 + 2P( 



Psb-P22Psb-P2QP,b 



2PI2PQQ - 2P22P20P02P00 - Poo ' 

' P22P20P,h - P20P02Psb ~ ^20 A)2Ab + 2PooP,b + 2P22PooP,h + P2oPooPsb 
' P22P20P00Psb + P20P02P00Psb ~ PooPs^b " ^22^-00^6 - + ^'22^sf 



- P02Psi - P22P02Psf - P20P02Psf + ^ 22^ 20^02nf ' 



- ^'22^'ooi'sf + P02P0QPsf + P22Po2PooPsi 



-^20^'02^sf- 



(5.43) 



It can be checked that this equation boils down to the expression ( |4.15| l in the single-step 
memory approximation. 

The validity of equations (5.17i, (5.29 1 and (5.38 1 has been checked by comparison 
with diffusion coefficients calculated from direct numerical simulations of the corresponding 
persistent random walk processes. 



5.4. Triangular lattice 



Consider finally the triangular lattice with 6-fold symmetry. We denote the relative directions 
by numbers from to 5, following the conventions shown in figure |2|d). 

The symmetry of M is similar to the previous subsections, so that the structure of the 
problem is by now clear Having identified the matrix J and invariant measure P in equation 
\5.A\, the velocity auto-correlation is found to be 
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2m 



A") 



'13,16 

An) 



An) 

n,9 



'13.23 



19.16 



'^19,23 ^ "^19,30 ^^'^19,31 
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Letting ([) = exp(2;r//6), we find 
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which is again equation (5.6 1. The expression of the diffusion coefficient is thus given by 
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([53. 



6. Two-Step Memory Approximation revisited 

As seen in the previous section, the symbolic computation of \5A) quickly becomes tricky. 
However, an alternative to the above scheme can be found, provided that the walk has special 
symmetries. Returning to \5.2\, we write 



X 



i^;o,S"'-''"-2t;o) 



(6.1) 



xP(S"^'2-i;o|S"-i;o,'yo)p('yo,S"-yo), 

where we introduce the compact notation S' = R'T, and sequences in the exponent denote 
multiple composition: S'' ■ = S'" o S'"-' o • • • o S'', where each /<. takes values between 1 
and z- Note that in general S-' o S' 7^ S'+ ' when T is non-trivial. The transition pro babilities 
f (S-'2?;o|S-"t>o,5-"'i>o) can be seen as matrix elements Qj2-/i.ii-;o' ^° '■h^'- '6.1 1 may be 
rewritten as 



{vq-v„)^Y. L ^o-S'' '"'i;oQi„,,„_, •••Q^2,,■l;'(^^o,S''^'o)• 



(6.2) 



We would like to rewrite this expression as a matrix product. However, this is in general 
not possible, and further approximations are needed. Thus, assuming that the scalar product 
v„ ■ V{) factorises as 



and defining 



k=\ 



Qij = Qi.jS-'vQ ■ vq = P(S-''''vo|S^i;o,'yo)S-'fo • vq , 



(6.3) 



(6.4) 



equation (6.2 1 becomes 



/t nn- 1 J 



(6.5) 



where we have introduced the vectors V,„ = vq ■ S'"Vo and p,| = p{vo,S'^ vq). As can be seen. 



equation ( 6.5 1 has an appropriate matrix form and can easily be resummed over n to compute 



the diffusion coefficient ( |2.1[ ). 

Since Q is a z x z matrix, equation (6.5 1 is much easier to evaluate than (5.4 1. The 
trouble is that equation ( |6.3| l is in general incorrect, and turns out to be strictly valid only 
for one-dimensional walks. Nonetheless, it may also be applied to higher-dimensional walks 
satisfying special symmetry conditions. We consider the different geometries separately in the 
following and discuss the conditions under which equation ( |6.5| l can be applied. For higher- 
dimensional lattices, we recover by this simpler method the results obtained earlier under the 
relevant symmetry assumptions. 



6.1. One-dimensional lattice 



The result (5.17 1 follows from equation (6.5 1. Indeed, R'"'- ''ivq • vq = R'"'' ^''vq • = ±1 
according to the parity of /„ + • • • + ii, and since this is also a property of the product 
n^=i R'^i'o ■ Vq, we see that equation (6.3 1 is valid. 
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The vector p, on the right-hand side of equation \6.5\ is 



Pi =p{-v^v) 

P2 = p(v,v) = 



2(l-Poo+Ao) ' 
Pw 



2(l-Poo + Ao) ■ 
The vector V,, on the other hand, has components 

Vi = -1, 

V2= 1. 



(6.6) 



(6.7) 



The matrix elements Q,- ^ are defined according to equation (6.4 1, 

Q = 



-Pn Pm 
-Pw Pqq 

Pio-1 1-Poo 
-Pu) Poo 



(6.8) 



Considering equation (2.1 1 and plugging the above expressions into equation (6.5 1, we 

have 

£'2sma = £'nma{i+4V+[I2-Q]"'p} , (6.9) 
and we recover equation ( 5. 17 1. 

6.2. Two-dimensional honeycomb lattice 



Consider equation (6.3 1 in the case of a honeycomb lattice. The operation S'v is a clockwise 
rotation of v by angle if i = 1, 7r/3 if i = 2, or tt if / = 3. The operation S'' '"v is thus 
a rotation by angle [2(/i + ... + !„) — 3n]7r/3, and the scalar product 

S'l v • V = cos[2(/i + ... + /„)- 3n]7i;/3 . (6.10) 

This expression is, however, in general different from the product 

S'"v ■v-S''v-v = ]~[cos[2(4) - 3]7r/3 . (6.1 1) 

k=l 



This is so, for instance, when n — 2 and ii = i2 = 1, for which (6.IO1 yields —1/2, whereas 
( [6TT] ) yields 1 /4. 

There is however a special case under which the product structure that we seek can be 
retrieved, as follows. There are a priori nine transition probabilities P{S-'''v\S-'v,v). There 
are, however, a number of left-right symmetries in the system which reduce the number of 
independent transition probabilities to three: 

P{S^-\\S^v,v), P{S^-\\S\,v), P{S^'^v\S^v,v). (6.12) 

In the event that the two probabilities f (S''^t>|S't>, t;) and f (S'''t>|S't>,i;) are equal, 

P{S^-^v\S^v,v) = P{S^-^v\S^v,v) = Pss, (6.13) 

which is to say that forward-left and right scatterings are treated as identical events, and the 
number of independent parameters reduces to two, which we take to be Poo and P^s- 
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In this case, the expression of the diffusion coefficient can be obtained in a way similar 
to equation (5.17 1 for the one-dimensional lattice. This is so because 



V - v = cos 



;;i...,i„_i.2 



V V = COS 



V - V — COS 



{[2(/i + ... + /„_i + l)-3n]|} 
COS - COS I [2(/i + . . . + /„_i) - 3(n - 1)] 3 j 
+ sin ^ sin | [2(/i + . . . + /„_i ) - 3(« - 1 )] || 
{[2(/i + ... + /„_i+2)-3«]|} 
|cos{[2(/i + ... + /„_i)-3(n-l)]|} 
in - sin I [2(/i + . . . + /„_i ) - 3(n - 1 )] - j 
{ [2(/i + ... + /„_!+ 3) -3«] I } 
{[2(/, + ... + /„_i)-3(«-l)]|}. 



(6.14) 



= —cos 



(6.15) 



(6.16) 



Thus, given the symmetry between forward-left and right scatterings, the two sine 

5 as the cc 



contributions in equations (6.14i and (6.15 1 cancel, whereas the cosines add up to 1: 

= S"-' 



'"-^v-v{S^v-v + S-^v-v), (6.17) 



"-^v ■ vS V ■ V . 



(6.18) 



We therefore retrieve an effective product structure, as in equation (6.3 1, and can compute the 
diffusion coefficient using equation (|6.5|l, with 



1 



and 



V: 



Q 



3(2-2A,-Poo) V l-2As 
1 

-1 



Pss I/2P0O-I/2 

-Poo 



(6.19) 
(6.20) 

(6.21) 



We obtain the expression of the diffusion coefficient for the symmetric [in the sense of 
equation ( |6.13| l] two-step memory approximation on the honeycomb lattice; 

3(1 -Poo) (2Pss+^'oo) 



D 



2SMA 



(3 -4ns + Poo) (2 -2ns -Poo) 



(6.22) 



This is equation (5.33 1. 



6.3. Two-dimensional square lattice 

For the two-dimensional square lattice, recall that T is the identity and the operation S'v = R'v 
is a anticlockwise rotation of v by angle i7z/2, with / = 0, . . . ,3. The operation S'' ■ '"v is thus 
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a rotation by angle {i\ 



i„)n/2. Equations similar to (6.14 1-(6.16 1 hold§| 



R'i-^'"vv= cos 



(/i + ••• + /„) 2 



cos — cos 
2 

— Sin — sin 
2 



(/i H h/„-i)- 

(i'l H l-/„-i): 



R' 



1. ..,(„_[ 



V ■ V R'" V ■ V 



(5,„,i -5,„3)sin (/iH h/„-i): 



The last term drops out provided 

Under the additional assumption that 

P{R^+'v\R^v,v) ^ P{R^+'v\R\,v) , 



(6.23) 
(6.24) 
(6.25) 



we retrieve an effective factorisation similar to that postulated in (6.3 i, and we can then use 
(6.5 I to obtain the corresponding diffusion coefficient. We refer to equations (6.24 1 and (6.25 1 
as defining a complete left-right symmetry. In this case, the invariant distribution is the 
solution of 



p{v, R'v) = Y,P{R'+jv\Rjv,v)p{R'v, Wv) , 



J^p{v,R-v)^-. 



(6.26) 

(6.27) 



We solve these equations for p{v,v) and p{v, —v), identifying p{v, R'v) and p{v, R^v), and 
define 



V = 



Q = 



with the transition matrix 

-P{v\-v,v) P{~v\v,v) \ / -P22 Pq2 
-P{~v\^v,v) Piv\v,v) J y -Pio Pqo 

The resulting expression of the diffusion coefficient is identical to (5.43 1. 

The validity of this expression extends to (/-dimensional orthogonal lattices under the 
symmetry assumptions ( |6.24[ )-( [6.25[ ). 



(6.28) 



(6.29) 



7. Conclusions 

We have shown that it is possible to find exact results for the diffusion coefficient of persistent 
random walks with two-step memory on regular lattices, by finding the matrix elements which 
give the velocity auto-correlation function and then resumming then. 

§ Note that, in general, we have the decomposition 

R'l ^,. I, = ^ sgn(tBi, . . . , ti)„)R'' u- ■ ■ ■ R'"i'- I'm,, , 
ffli (Bi,e{o.l} 

where we introduced the notation v^i = v if CO = and Vm = Vj_ if 0) = 1, and the function sgn(Q)i , . . . , (0„) = ±1, 
depending on the sequence (Bi , . . . , (B„ . Equation l |6.5^ would then be replaced by a more complicated expression 
involving the mixed products of two matrices P{R^'^'v\RJv, u)R^d ■ v and P{RJ'^'v\R^v, v)RJv ■ v±. 
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We have applied the resuhs obtained here to approximate the diffusion coefficients of 
certain periodic bilUard tables in II H . 

The extension to lattice random walks with longer memory is possible, albeit difficult for 
obvious technical reasons. Finally, we remark that the extension to lattices in three dimensions 
is not direct, since in that case, additional information must be specified in order to uniquely 
define relative directions 1 12]. 
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Appendix A. 2SMA on a honeycomb lattice 



In analogy to equation (5.13 1, we may write 
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We have the three sets of equations 
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Proceeding with our analogy, we seek linear comb 
matrices on the left-hand side of the above equations, so as to obtain a single matrix equation 
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nations of the elements in the rows of the 
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similar to equation ( |5.15| l. Considering the elements in equation ( |A.2[ ), we write 
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(A.5) 



Comparing with equation (A. 2 1, we infer 

[ = exp(-2;;r/3) = -^-^ 
Applying the same procedure to equations (|A.3[)-(|A.4|i, we obtain 



= exp(2;7r/3) = -^-^ , 



(A.6) 




(A.7) 



The system of equations ( |A.5| l, ( |A.7| i 
equation 



(A.8) 



|A.8|l reduces to the single recursive matrix 
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^00 

fPoi f-Pn (p^Pii 
Appendix B. 2SMA on a square lattice 




(A.9) 



In analogy to equations ( |A.2| i-( [A~4| ), there are 64 different entries among the 256 elements of 
M", which can be obtained through the set of equations 



/ a 
d 
c 



(" 

d<" 
"3 

^3 

V bi" 

(n 



d^" 

^7 

a'" 

^8 
d^" 

Jn 



d« 



°9 



-10 



b^" 

°10 



(") 

16 
(«) 
16 
in) 
16 
(«) 
16 



°13 

An) 
^13 

Dj3 



°14 

a 14 

^(«) 

°14 



°15 

dW 
°15 

An) 
'-15 



\ 



°15 



Combining these quantities, we let 
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provided 



^^ = l^^,, = exp{ik7!:/2), A; = 0,1,2,3, 
We have the following identities, 
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